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STUDIES 


THE  UTILIZATION  OF  THE  ANALYTIC  GBAPH 

METHOD  TO  DETERMINE  THE  BEST  SHAPE  OF 

f-L0VJ 

THIN  DELTA  WINGS  IN  SUPERSONIC  -GtfftRENTD* 


by  ADRIANA  NASTASE** 

C.  NASTASE*** 

FRANC INE  ATANASIU- MOLDOVAN**** 


In  the  present  work  an  analytic  graph  method  is  used 
( bes.tJ 

(by  the  first  author)  t^determlne  the  best  shape  of  a thin 
delta  wing  from  among  the  original  thin  delta  wing  class  In 
which  the  distribution  of  the  vertical  speed  of  disturbance 
Is  expressed  In  the  form  of  second  and  third  order  homogenous 
polynomials.  In  addition,  the  thin  delta  wing  satisfies  and 
meets  the  geometrical  and  aerodlnamlc  conditions  of  natural 
law. 


1.  INTRODUCTION 

The  analytic  graph  method  Is  utilized  In  the  present  work 
(1)  for  the  effective  determination,  at  a given  Mach  cruising 
speed  M,  the  best  optimum  form  of  thin  delta  wings  with  sub- 
sonlff'rBoges , which  belong  to  the  original  class  of  thin  delta 

1. 


wings  and  which,  in  addition,  satisfy  the  following  conditions 
of  law:  the  lift  power  and  moment  of  pitch  are  given,  and  the 
- . Wf* ** *** ****^disturbance  speed  u is  finite  aLong  the  subsoni^'e^ge^ 
of  the  wings,  in  order  to  avoid  the  formation  and  falling  of 
the  vortexes  which  have  the  tendency  to  appear  along  that  edge 
at  the  Mach  number  of  cruising  speed  M.  The  incidents  of  the 
thin  delta  wing  are  presupposed  to  be  sufficiently  small,  in 


such  a manner  that  we  are  able  to  apply  the  results  of  disturbance 
theories . 

The  vertical  speeds  of  disturbance  on  the  wings^Wjls  analyzed 
by  using  second  and  third  order  homogenous  polynomial  tests* 


* The  authors  thank  Mrs.  Denise  Vallee-Gulraud , main 
scientific  researcher  and  Mr.  M.  Bols,  scientific  researcher 
at  ONERA  (France)  for  valuable  advice  in  laying  out  the  calculus 
program. 
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Thus  we  can  anply  the  results  of  P.  Germain’s  higher  order 
cone  flow  theories,  (3). 


The  axis  of  disturbance  speeds  u on  the  wing  Is  obtained 
with  the  help  of  the  hydrodynamic  analogy  method,  through  the 
suoerimDOSltion  of  the  contributions  of  the  borders  and  edges 


on  the  wing  (4),  (5)* 


The  study  of  the  optimum  form  of  a thin  delta  wing  through 

the  analytic  graph  method  leads  to  the  determination  of  the 
V 

optimum  values  /0  of  the  oarameter  of  similitude  / = B1  of  flow, 

v 

which  minimizes  the  wave  resistance  expression,  = f(^)»  For 
that  purpose  an  Imaginary  transformed  thin  delta  wing  (projected 
In  a fixed  plane)  is  used,  which  Is  the  only  function  of  the 


parameter  of  similitude  /. 


The  transformed  thin  delta  wing  is  placed  In  a suitably 
chosen  supersonic  flow. 


2.  Axial  and  Vertical  Disturbance  Velocity  Expressions 


The  thin  delta  wing  is  referred  to  a trlorthogonal  system 
of  axis  Ox^gX^,  which  has  the  origin  In  the  peak  of  the  wing, 
axis  Oxj  parallel  with  the  velocity  at  Infinity  UU*  again  the 
wing  Is  plotted  and  numbered  in  plane  (fig»l)» 


The  vertical  disturbance  velocity  w„  =- 


on  the  Initial 


thin  delta  wing  Is  presupposed  as  expressed  In  the  form  of 


, m.naiWrv 


second  and  third  order  homogenous  polynomials 

*.  “ -jj-  * *»(«>io  + Woilyl)  + + «®nlyl  + M>oa 3/‘).  (1) 

In  the  following  similarity  transformations  (1), 


plane  0x^21  the  designated  transformed  plane , corresponds 

to  the  initial  nlain  Ox, x . 

1 2 

The  initial  thin  delta  wing  projection  in  plane  Ox^x2 
is  an  Isosceles  triangle  0A^A2  in  which  the  height  corresponding 

l . 

to  peak  0 is  h^,  and  its  base  is  21|(fig.l). 

In  the  succeeding  similarity  transformation  (2)  the 
transformed  thin  delta  wing  projection  in  plane  ^x^x^  1®  an 
isosceles  triangle  "^A^X?  which  has  height  ! and  base  2 (fig. 2), 
and  the  vertical  disturbance  velocity  wg  = vfa  acquires  the 
following  form 

*’•  * fi'  ~ + ®oil  y I)  + + u>„  |y  I + «>0ly*).  (3) 

Coefficients  Wjj  andw'jj  from  (!)  and  (3)  are  corinected 
through  the  relations: 


7 

Fig.  1 


The  expression  of  the  disturbance  axis  velocity  u on  the 
initial  thin  delta  wing  with  subsonic  leading  edges,  corresponding 
to  the  distribution  (1)  of  the  vertical  disturbance  speed  w, 
is  of  the  form: 


Fig.  2 


In  a similar  manner,  the  expression  of  the  disturbance  axis 


The  constants  A^j  and  Cjj  of  the  axis  of  disturbance 
velocity  on  the  transformed  thin  delta  wing  are  connected  by  t 

coefficients of  the  vertical  disturbance  speed  K on  that 
wing  through  homogenous  and  linear  relations  deduced  from  the 

I : r 

compatibility  conditions  of  P,  Germain  (3): 


In  = 58  + «•»’  w.i. 

(io», 

^M  = a|i'©„  + «?,'»•„ 

(10b| 

= a0e  + •fi'Su  + *8  *®M> 

(10ci 

In  = 58  + 5j?»u  + 58  u>„, 

(lOd) 

Z/M  — e{l  m>  u. 

(lOe) 

6. 

i 


v>  l 

If  y = BJ  as  noted  earlier  and  the  eliptlc  Integrals  of 
the  first  and  second  cases  E(k)  and  K(k),  by  the  mode 

k = H~-~A  (11* 

then  the  coefficients  «i}’t  which  Intervene  in  the 

expressions  (10  a,  b,  c,  d,  e)  of  the  constants  of  the  axis  of 
disturbance  velocity  u on  the  transformed  thin  delta  wing  with 
subsonic  leading  edges  are  , respectively,  of  the  form  (6),  (?),. 
(8): 


54  - • 

-«» 


m 


_ _ 2[E(k)  - v»JT(fc)] 


1- V* 


(12u 


(12b 


(12c) 


■S'  = 


2v*[E(k)  - K(k)] 
nN. 


(12d) 


The  following  expression  is  noted  with  Ngt 

^i  = (l-  2v*)  E(k)  + v*iC(fc),  ( 13) 

= Y [2(3  - 5v*  + v«)  E(k)  - v»(3-5v»)J[(i)].  (Ha) 

^ [<4  + V*)  E*(k)  - 8 v*J?(fc)  (*)  + 3v«Jf*(fc)],  (14b) 

Y t(1  + V,)  E(k)  ~ 2 ( 140 
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(14d) 


offl  = - ~ [(4  - 7 v*  + v4)  E(k)  - 2(1  - 2v»)  *(*)], 

aii’  = — [(12  - 17 v*  + 10  v4)  E*(k)  - 

nN9 

- 4 v*(2  - v*  + v4)  E(k)K(k)  + v4(l  + 2v»)  *»(*)],  (14e) 

Sff  = “‘[2(2  - v*)J m -(3  -v*)  K(k)],  (14f) 

«8 ' = (14g) 

7C 

* - - . _ 

The  expression  noted  with 

A’,  = (4  — 19v*  + 4 v4)JB*(Jfc)  + 8v*(l  + »*)E(k)K(k)  - 5v4J£*(*).  (15) 

The  explicit  forms  of  the  relatlohs  of  connections  which 
occur  In  problems  of  the  optimum  transformed  thin  delta  wing 
are  the  following: 

- Lift  conditions  of  the  transformed  thin  delta  wing 
to  be  given  as 

C,  — 8 ^ Muljdxjdp  --  AI0wlu  + A21«>0,  A.mwm  + 

-\  A31m)u  + A33woa  = 2*  [(2 alSi  + «8l  )®jo  + 

l) 

+ (25^'  + «i¥  )w».]  + [(23U)  + SIS  )w„  + ^20oV  + 

+ H — g'-j  + (2flo*’  + O**1  (M’ojj  — Vi%  = — - ■ O'  ( / & 3 

The  moment  of  pitch  condition  of  the  transformed  thin 
delta  wing  to  be  given  as: 


8 


ft. -8  j«a3jfd21«iy=r,uwJ0+  f3lw01  4 TMW.M  +rillw1,  + 


I 

V 

\ 

i 


3 _ -w  4 ^ ^ 

•(-  I 3a:«Vi  — (-^ao^io  4"  Aj,t0#1)  -f  ^ (A30u>M  + A3,tt5,i  4~ 


~^32wm)  — V„t  = - ' • 


07  ) 


The  disturbance  axis  velocity  condition^  is  to  be  finite 
alone  the  subsonic  leading  edgeSof  the  transformed  wings  (1) 
leading  in  this  case  to  the  relation  of  connection  of  the  form: 

/’jia  j420  + -4  j2  — 0, 

jfj  — ^30  4"  ^32  — 0. 


. . (iraS) 

(I*- 

(li'1  f / 9 >6  3 


If  relations  ( 10a,b,c,d ,e)  are  taken  into  consideration, 
these  relations  are  written  in  the  form: 

7 — ^20*^10  4*  4 21^01  — ( <*00  4"  <*20  )W>J0  4"  (Oi®1  4 ®M  )tCgl  ~ ®>  H'* 

^ 3 — ]-  T3iwlt  4-  T3Iu>qi  = ( «oo  4-  oiJ1  )wt0  4 

4-  (<*oY  4-  aiV  )«ji  4-  (215  -f  25  )2>0I  = 0.  (-' 

The  expression  of  the  coefficient  of  wave  resistance  of 
transformed  thin  delta  wing  is  in  the  form: 

5«5A  dJi«ly  - “*«•«&  4-  (fi„l04- 

+ Smoi)»i«»»m  4-  Q||||4  4-  <ntJM  4-  nil00)^ltw01  4- 
4-  (fl*3»i  4-  Uitio)v>9lwto  4-  ( ft,,,,  4-  Ww„)»#1i5a  4. 

4 (2M1,  4-  Qnit)trat»w  4-  UjMoi&  4 nai]tof,  4 

9. 


<7  70 

the 


f 


i 


1 


+ a33M«4  + (liJ3io  + i^aoi)M,aoS11  +(033ao  + 
+ i^330»)wmwua  -f  (ii3321  f iiyJli)wnWUi, 


(»-!)■ 


(21) 


Continuing,  the  following  definite  Integral  Is  marked  with  3k 


3k 


i 


3U  — 


yk&y  t 

u -r 

_n(2<) ! 
2a  + ,(t!)s 


n 2 «« !)* 
(2t  + l)! 


(22a) 

(*  = 21),  (22b) 
(it  = 2<  + 1).  (22c) 


From 

which 

for  k = 0 

,1,2, 3, 4, 5 

the  following  values  for  g 

result: 

n _ * 

3o-~> 

3x  =1, 

(23a) 

3 ,7t 

2,  = 

3 

(23b) 

3«  = — . 

16 

3,=  -• 

15 

(23c) 

With 

these 

notations 

linmkj,  the 

constants,  which  appear  In 

exoresslon  (21)  of  wave  resistance^  of  the  transformed  thin 


delta  wing  can  be  written  under  the  form 

Qtmkl  = — (<i't)'3k  + a'i>3klt) 

to  + 2 


(24 


(to  = 2,3,  it  = 0,1,...,  (to  - 1),  j =0,1  ) 


and,  respectively, 


°'mki  ~ mTa  [<>Sk  + + irfy ) Sk+i] 


(AS) 


(w  = 2,3  It  = 0,1,. j = 0,1,2). 


To  explain  by  means  of  example,  in  this  case,  the  coefficients 


are  in  the  form: 


^imo  — ~ (2a(K)  + <*!»),  ■ ^mio  — ~ (3flo#  + 2ag,’),  (26u 


nK01  - S (2aiV  + 3g),  Oim  = -I  <33fi’  + 25g>)  (261 


for  n = 2 and  m = 2 , 

a 2tc 


e”  (^aoo  + aw)>  ^uio  — + 2a^)y  (27a  t 

u 15 

= 2f  (2S«  + 52’),  nM11  = 1 (35g  + 2ag'),  (271, 

®*uo  =»  — (43$  -f  3 ajSl),  (27c) 


^uu  — — (4aJ,V  -f  3ag') 

for  n = 2 and  m = 3, 

O — 

“,*o.  = — (2S«i  + aft), 


S.»0  = A (3«iS*  + 2aS),  (28a) 


fi„oi  = ^ « + 331V  + okV),  S,«u  = ^ K63g>  + 45ff  + «). 


^ a/A^rifl 


Bum  = 2k-  (25&>  + 31V),  ft,*.  - — (32fi  + 22«,)  * <28c) 

for  n = 3 and  m * 2 respectively, 

«moo  - V <25&  + 5">»  “»»•  “ “ (3S&  + 2«»)>  C2®*) 

5 9 

&,»,  = (23)3  + 3ffi),  nM11  = A (631V  + 4«ff  + eg'),  (29b) 

o 9 

0MM  = A <45&  + 3S&),  (2»d) 


11 


r 


5*»*i  “ ~ (2«ai*'  -f  15«a»  + 3cfi'), 
60 


5**»*  — “ (4ao*  + 30^’ 

1J 


(29e) 

(29f) 


for  n = 3 and  m = 3* 


The  equations  (1),  (9)  obtained  through  canceling  all 

L *SkJ 


( 

of  the  variation  of  coef  f 1 cl  entlf^qRre  ^in  thl 
form: 

2ilril)0wl0  4-  (H  tsoi  + 5**1*)  «>n  + (5,,#,  + ft 


s case  of  the 
> 


»)  + 


+ (H***l  + ^Ull)  WII  + (11**0*  4-  5****)  »**  + 4- 

4-  X«*'PM  4-  Xtflo=0,  (30) 

(H**l*  + H**0l)  *®M  + 2niin«014-(Q  Ml*  + H**l*)  *>**  + (11**11  + 

4-  H**n)  wn4"(ll*»i»  f U**ii)  w*»  + 4-  X1*’!''**  4- 

4-  X,  T„  = 0, 


(31) 


(O2M0  4-  Qjioo)  »io  4-  (11**01  4-  Umi*)  wn  4-  25****  *>u  4* 

4-  (11**01  4-  H3,10)  «>n  4-  (11*»**  4-  11****)  •*’**  4-  ^<nA*o  4- 
4-  X^r*,  4-  X*Pm=0,  (32) 

(H*3]o  4-  11**01 ) wi*  4-  (ll**n  4-  5,1U)  P’01 14-|  (5U„  4-  5***j)  »*•  4* 
4-  25WI«5„  4-  (Omi,  4-  nMM)  5f„  4-  X'»A„  4-  X«**fM  4- 

s m A 


4-  X.T,,  = 0, 


(33) 


(ll***o  4-  n*MI)  «?i*  4-  (Q**u  4-  5„„)  icn  4-  (5*»» * 4-  5mm)  «cm  4- 
4-  (5UI1  4-  5,m)  «n  4-  2Q****iOoi  4-  X(,IAM  4-  4- 

4-  \f  „ = 0. 


(34) 


The  equations  (30  - 344,,  together  with  the  relations  of 
connection  (16),  (17),  (I8a,b),  (19)  and  (20),  formH^ly^fem 
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I 


of  algebraic  equatl  ms  from  which  Is  determined  optimum  values 
of  the  coefficients 'wjq.w  /w20,wxl  and  ff02  of  the  vertical 

velocity  of  disturbance  w , as  well  as  the  values  of  the 

and 

multipliers  X'11,  X<*’,  x,  \|i,  x»  f or  a given  value  of  the  parameter  of 

v / 

similitude  / = BJi. 


2ft»to  ^uci  ^Jttt 

K 

C 

1 

0 

^ nu  Qooo 

2ft w, t ft  lit  ftono 

ft oi,  ft 30,0 

ft  1300  ft 3101 

A„ 

n, 

0 

0 

n • n 

HOt  3000 

W*t  20m 

niin3W 

ftootoftjuo 

rM 

0 

t 

n • 0 

nw  sot 

ftoMAo* 

20 

JJU 

n • o 

JMt  JJ/f 

\ 

C 

0 

0 

n - n 

0300  3100 

203301 

K 

G 

0 

0 

\ 

\ \ 

K 

K 

0 

0 

0 

0 

r. 

r„  rM 

r , 

r„ 

0 

0 

0 

0 

L 

Tn  0 

0 

0 

0 

0 

0 

0 

0 

o T„ 

T* 

T» 

0 

0 

0 

0 

APPENDIX  1 


any 

If  the  order  of  unknowns  X<»,  >. 

Is  chosen  presupposing  the  following  order  of  equations  In 
systems  (30)  - (34),  (l6),(17),  (19)  and  (20),  then  the 
quadratic  matrix  for  rank  m = 9 of  the  system  of  coefficients 
formed  from  the  linlar  equations  mentioned  above  has  the  form 
In  appendlx(l),  and  different  Independent  single  terms  for  zero 
appear  In  the  sixth! 1 6 ) and  seventh  (1?)  equations  and  have 
the  values  respectively.  Because  the  solution 

V V 

of  the  system  is  nc* changed  for  diverse  values  of  the  Initial 
givens  (Cj^  and  Cm  ) the  matrix  of  Independent  terms  Is  considered 
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to  be  a rectangular  matrix  R of  the  tyoe  m x n = 9 x 2, 


f -®ll 

J2U' 

1 R* 

Ru 

Rn. 

(3:.,  ( 3 


In  which  R = 0 with  the  exception  of  the  tefms  Rfi  = 1 
and  R72  = 1.  If  the  solutions  of  the  system  obtained  for  the 
two  sets  of  values  for  liniar  terms  R^  and  R i is  noted  with 

tf 

wij  an<^  wij  ' * = * * then  the  optimum  values  of  the 

coefficients  w^j  have  the  values  given  by  the  relation 

w*l — w*i  + v~  • (3  ( J6,  ) 


For  the  determination  of  the  optimum  form  of  projection  in  a 
plane  at  the  Mach  flight  number  of  M =2,  the  range  of  admissible 


values  of  the  span  of  the  wing 


/ ^21 

Lx  = - 


h.. 


the  range  of  values 


contained  between  0,2<{<o,55  Is  considered  ( the  value 

$ = 0*5774  corresponds  to  the  sonic  leading  edge  for  Mach  cruizing 


speed  M = 2). 


Therefore,  the  following  range  of  values  admissible  for 
result  : 

0,3464  < v < 0,  9626.  (37) 

V 

A set  of  values  of  the  parameter  of  similitude  y contained  in 
this  range  of  admissible  values  is  considered,  with  the  help  of 


Bfifu 


itlHfflimilg 


14* 


T 

If&rhelght  Is  calculated) 

the  following,  r2i)Vand  the  curve  (C^7 eight  = f(v)  (for 

opr 

Q = constant)  which  has  the  asoect  in  figure  3,1s  drawn. 


With  the  help  of  this  curve  the  optimal  value  of  the 
parameter  of  similitude  y is  graphically  determined  (for  which 
this  curve  touches  its  minimum). 


\>  V 

The  optimum  value  jfQ  of  the  parameter  of  similitude  y, 

together  with  the  given  area  of  projection  In  the  plane, 
determines  the  optimum  form  of  the  initial  thin  delta  wings' 
projection  on  the  plane.  ♦ 


r* 

Hm‘* 

1* 

iV* 

I 3 

i 

! 

i 

i 

1 

1 

1 

1 

1 

i'll  1 

! 

i 

1 

1 

1 

i 

1 

i ] 

i 1 

i 

i i!  i 

— wri  1 

Vqmi  qs 

a*  i 4* 

Liilk  * 
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To  make  it  clear,  consider  if  you  will,  the  following 
given  initially: 

- the  number  of  Mach  speed,  * 2, 

- the  coefficient  of  lift,  = 0,3, 

- the  coefficient  of  the  moment  of  pitch,  Cm  = 0,23, 

- the  area  of  projection  the  plane,  A0  * 75  m2, 


ft 


the  optimum  form  of  the  projection  on  the  plane  of  the 
initial  thin  delta  wing  Is  the  isosceles  triangle  which  has 
a height  of  h1  = 12,63  m and  a base  of  = 11,86  m,  and 

Mm 

the  optimum  unit  span  is  M- - » 0.47. 

hi 

The  minimum  value  of  wave  resistance  of  the  initial  thin 
delta  wing  is  given  by  the  relations 

Ct  = icd  = 0,327. 

* 

Accepted  by  the  editors  on  Dec.  29,  1972 
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